Optimizing the frame duration for data-driven rigid motion estimation in brain PET imaging.

PURPOSE: Data-driven rigid motion estimation for PET brain imaging is usually performed using data frames sampled at low temporal resolution to reduce the overall computation time and to provide adequate signal-to-noise ratio in the frames. In recent work it has been demonstrated that list-mode reconstructions of ultrashort frames are sufficient for motion estimation and can be performed very quickly. In this work we take the approach of using image-based registration of reconstructions of very short frames for data-driven motion estimation, and optimize a number of reconstruction and registration parameters (frame duration, MLEM iterations, image pixel size, post-smoothing filter, reference image creation, and registration metric) to ensure accurate registrations while maximizing temporal resolution and minimizing total computation time. METHODS: Data from 18 F-fluorodeoxyglucose (FDG) and 18 F-florbetaben (FBB) tracer studies with varying count rates are analyzed, for PET/MR and PET/CT scanners. For framed reconstructions using various parameter combinations interframe motion is simulated and image-based registrations are performed to estimate that motion. RESULTS: For FDG and FBB tracers using 4 × 105 true and scattered coincidence events per frame ensures that 95% of the registrations will be accurate to within 1 mm of the ground truth. This corresponds to a frame duration of 0.5-1 sec for typical clinical PET activity levels. Using four MLEM iterations with no subsets, a transaxial pixel size of 4 mm, a post-smoothing filter with 4-6 mm full width at half maximum, and averaging two or more frames to create the reference image provides an optimal set of parameters to produce accurate registrations while keeping the reconstruction and processing time low. CONCLUSIONS: It is shown that very short frames (≤1 sec) can be used to provide accurate and quick data-driven rigid motion estimates for use in an event-by-event motion corrected reconstruction.

Ultra-Fast List-Mode Reconstruction of Short PET Frames and Example Applications.

Standard clinical reconstructions usually require several minutes to complete, and this time is mostly independent of the duration of the data being reconstructed. Applications such as data-driven motion estimation, which require many short frames over the duration of the scan, become unfeasible with such long reconstruction times. In this work, we present an infrastructure whereby ultra-fast list-mode reconstructions of very short frames (≤1 s) are performed. With this infrastructure, it is possible to have a dynamic series of frames that can be used for various applications, such as data-driven motion estimation, whole-body surveys, quick reconstructions of gated data to select the optimal gate for a given attenuation map, and, if the infrastructure runs simultaneously with the scan, real-time display of the reconstructed data during the scan and automated alerts for patient motion. Methods: A fast ray-tracing time-of-flight projector was implemented and parallelized. The reconstruction parameters were optimized to allow for fast performance: only a few iterations are performed, without point-spread-function modeling, and scatter correction is not used. The resulting reconstructions are thus not quantitative but are acceptable for motion estimation and visualization purposes. Data-driven motion can be estimated using image registration, with the resultant motion data being used in a fully motion-corrected list-mode reconstruction. Results: The infrastructure provided images that can be used for visualization and gating purposes and for motion estimation using image registration. Several case studies are presented, including data-driven motion estimation and correction for brain studies, abdominal studies in which respiratory and cardiac motion is visible, and a whole-body survey. Conclusion: The presented infrastructure provides the capability to quickly create a series of very short frames for PET data that can be used in a variety of applications.

Real-Time Gain Control of PET Detectors and Evaluation With Challenging Radionuclides.

Accurate gain control of PET detectors is a prerequisite for quantitative accuracy. A shift in the 511 keV peak position can lead to errors in scatter correction, degrading quantitation. The PET detectors in a PET/MR scanner are subject to thermal transients due to eddy currents induced during gradient-intensive MRI sequences. Since the gain of silicon photomultiplier-based detectors changes with temperature, good gain control is particularly challenging. In this paper we describe a method that utilizes information from the entire singles spectrum to create a real-time gain control method that maintains gain of PET detectors stable within approximately ±0.5% (±2.5 keV) with varying levels of scatter and in the presence of significant thermal transients. We describe the methods used to combine information about multiple peaks and how this algorithm is implemented in a way that permits real-time processing on a field-programmable gate array. Simulations demonstrate rapid response time and stability. A method ("virtual scatter filter") is also described that extracts unscattered photopeak events from phantom data and demonstrates the accuracy of the photopeak for various radionuclides that emit energies in addition to the pure 511 keV annihilation peak. Radionuclides 52 Mn, 55 Co, 64 Cu, 89 Zr, 90 Y, and 124 I are included in the study for their various forms of spectral contamination.

PET Image Quality Improvement for Simultaneous PET/MRI with a Lightweight MRI Surface Coil.

Background During simultaneous PET/MRI, flexible MRI surface coils that lay on the patient are often omitted from PET attenuation correction processing, leading to quantification bias in PET images. Purpose To identify potential PET image quality improvement by using a recently developed lightweight MRI coil technology for the anterior array (AA) surface coil in both a phantom and in vivo study. Materials and Methods A phantom study and a prospective in vivo study were performed with a PET/CT scanner under three conditions: (a) no MRI surface coil (standard of reference), (b) traditional AA coil, and (c) lightweight AA coil. AA coils were not used in attenuation correction processing to emulate clinical PET/MRI. For the phantom study, PET images were reconstructed with and without time of flight (TOF) to assess quantification accuracy and uniformity. The in vivo study consisted of 10 participants (mean age, 66 years ± 10 [standard deviation]; six men) referred for a PET/CT oncologic examination who had undergone imaging between October 2019 and February 2020. Assessment of image quantification bias (defined as the standard error of the mean values) was conducted by comparing mean liver region of interest standardized uptake values with the no-coil standard of reference. A Wilcoxon signed-rank test was used to establish significance. Results For TOF and non-TOF, respectively, the phantom study revealed a mean PET quantification bias of -9.0% and -8.6% with the traditional AA coil and a mean PET quantification bias of -4.3% and -4.0% with the lightweight AA coil. The coefficients of variation reduced from 4.3% and 6.2% with the traditional AA coil to 2.1% and 2.7% with the lightweight AA coil, which demonstrated a homogeneity benefit from the lightweight coil that was greater with, versus without, TOF reconstruction. For the in vivo study, the mean liver standardized uptake value error was -5.9% with the traditional AA coil (P = .002 vs no coil) and -2.4% with the lightweight AA coil (P = .004 vs no coil). Conclusion The lightweight anterior array coil reduced PET image quantification bias by more than 50% compared with the traditional coil. Using the lightweight coil and performing time of flight-based reconstruction each reduced the variation of error. © RSNA, 2020 Online supplemental material is available for this article.

Quantitative and Qualitative Improvement of Low-Count [68Ga]Citrate and [90Y]Microspheres PET Image Reconstructions Using Block Sequential Regularized Expectation Maximization Algorithm.

PURPOSE: There are several important positron emission tomography (PET) imaging scenarios that require imaging with very low photon statistics, for which both quantitative accuracy and visual quality should not be neglected. For example, PET imaging with the low photon statistics is closely related to active efforts to significantly reduce radiation exposure from radiopharmaceuticals. We investigated two examples of low-count PET imaging: (a) imaging [90Y]microsphere radioembolization that suffers the very small positron emission fraction of Y-90's decay processes, and (b) cancer imaging with [68Ga]citrate with uptake time of 3-4 half-lives, necessary for visualizing tumors. In particular, we investigated a type of penalized likelihood reconstruction algorithm, block sequential regularized expectation maximization (BSREM), for improving both image quality and quantitative accuracy of these low-count PET imaging cases. PROCEDURES: The NEMA/IEC Body phantom filled with aqueous solution of Y-90 or Ga-68 was scanned to mimic the low-count scenarios of corresponding patient data acquisitions on a time-of-flight (TOF) PET/magnetic resonance imaging system. Contrast recovery, background variation, and signal-to-noise ratio were evaluated in different sets of count densities using both conventional TOF ordered subset expectation (TOF-OSEM) and TOF-BSREM algorithms. The regularization parameter, beta, in BSREM that controls the tradeoff between image noise and resolution was evaluated to find a value for improved confidence in image interpretation. Visual quality assessment of the images obtained from patients administered with [68Ga]citrate (n = 6) was performed. We also made preliminary visual image quality assessment for one patient with [90Y]microspheres. In Y-90 imaging, the effect of 511-keV energy window selection for minimizing the number of random events was also evaluated. RESULTS: Quantitatively, phantom images reconstructed with TOF-BSREM showed improved contrast recovery, background variation, and signal-to-noise ratio values over images reconstructed with TOF-OSEM. Both phantom and patient studies of delayed imaging of [68Ga]citrate show that TOF-BSREM with beta = 500 gives the best tradeoff between image noise and image resolution based on visual assessment by the readers. The NEMA-IQ phantom study with [90Y]microspheres shows that the narrow energy window (460-562 keV) recovers activity concentrations in small spheres better than the regular energy window (425-650 keV) with the beta value of 2000 using the TOF-BSREM algorithm. For the images obtained from patients with [68Ga]citrate using TOF-BSREM with beta = 500, the visual analogue scale (VAS) was improved by 17 % and the Likert score was increased by 1 point on average, both in comparison to corresponding scores for images reconstructed using TOF-OSEM. CONCLUSION: Our investigation shows that the TOF-BSREM algorithm improves the image quality and quantitative accuracy in low-count PET imaging scenarios. However, the beta value in this algorithm needed to be adjusted for each radiopharmaceutical and counting statistics at the time of scans.

NEMA NU 2-2007 performance characteristics of GE Signa integrated PET/MR for different PET isotopes.

BACKGROUND: Fully integrated PET/MR systems are being used frequently in clinical research and routine. National Electrical Manufacturers Association (NEMA) characterization of these systems is generally done with 18F which is clinically the most relevant PET isotope. However, other PET isotopes, such as 68Ga and 90Y, are gaining clinical importance as they are of specific interest for oncological applications and for follow-up of 90Y-based radionuclide therapy. These isotopes have a complex decay scheme with a variety of prompt gammas in coincidence. 68Ga and 90Y have higher positron energy and, because of the larger positron range, there may be interference with the magnetic field of the MR compared to 18F. Therefore, it is relevant to determine the performance of PET/MR for these clinically relevant and commercially available isotopes. METHODS: NEMA NU 2-2007 performance measurements were performed for characterizing the spatial resolution, sensitivity, image quality, and the accuracy of attenuation and scatter corrections for 18F, 68Ga, and 90Y. Scatter fraction and noise equivalent count rate (NECR) tests were performed using 18F and 68Ga. All phantom data were acquired on the GE Signa integrated PET/MR system, installed in UZ Leuven, Belgium. RESULTS: 18F, 68Ga, and 90Y NEMA performance tests resulted in substantially different system characteristics. In comparison with 18F, the spatial resolution is about 1 mm larger in the axial direction for 68Ga and no significative effect was found for 90Y. The impact of this lower resolution is also visible in the recovery coefficients of the smallest spheres of 68Ga in image quality measurements, where clearly lower values are obtained. For 90Y, the low number of counts leads to a large variability in the image quality measurements. The primary factor for the sensitivity change is the scale factor related to the positron emission fraction. There is also an impact on the peak NECR, which is lower for 68Ga than for 18F and appears at higher activities. CONCLUSIONS: The system performance of GE Signa integrated PET/MR was substantially different, in terms of NEMA spatial resolution, image quality, and NECR for 68Ga and 90Y compared to 18F. But these differences are compensated by the PET/MR scanner technologies and reconstructions methods.

Clinical Evaluation of 68Ga-PSMA-II and 68Ga-RM2 PET Images Reconstructed With an Improved Scatter Correction Algorithm.

OBJECTIVE: Gallium-68-labeled radiopharmaceuticals pose a challenge for scatter estimation because their targeted nature can produce high contrast in these regions of the kidneys and bladder. Even small errors in the scatter estimate can result in washout artifacts. Administration of diuretics can reduce these artifacts, but they may result in adverse events. Here, we investigated the ability of algorithmic modifications to mitigate washout artifacts and eliminate the need for diuretics or other interventions. MATERIALS AND METHODS: The model-based scatter algorithm was modified to account for PET/MRI scanner geometry and challenges of non-FDG tracers. Fifty-three clinical 68Ga-RM2 and 68Ga-PSMA-11 whole-body images were reconstructed using the baseline scatter algorithm. For comparison, reconstruction was also processed with modified sampling in the single-scatter estimation and with an offset in the scatter tail-scaling process. None of the patients received furosemide to attempt to decrease the accumulation of radiopharmaceuticals in the bladder. The images were scored independently by three blinded reviewers using the 5-point Likert scale. RESULTS: The scatter algorithm improvements significantly decreased or completely eliminated the washout artifacts. When comparing the baseline and most improved algorithm, the image quality increased and image artifacts were reduced for both 68Ga-RM2 and for 68Ga-PSMA-11 in the kidneys and bladder regions. CONCLUSION: Image reconstruction with the improved scatter correction algorithm mitigated washout artifacts and recovered diagnostic image quality in 68Ga PET, indicating that the use of diuretics may be avoided.

PET Imaging Stability Measurements During Simultaneous Pulsing of Aggressive MR Sequences on the SIGNA PET/MR System.

The recent introduction of simultaneous whole-body PET/MR scanners has enabled new research taking advantage of the complementary information obtainable with PET and MRI. One such application is kinetic modeling, which requires high levels of PET quantitative stability. To accomplish the required PET stability levels, the PET subsystem must be sufficiently isolated from the effects of MR activity. Performance measurements have previously been published, demonstrating sufficient PET stability in the presence of MR pulsing for typical clinical use; however, PET stability during radiofrequency (RF)-intensive and gradient-intensive sequences has not previously been evaluated for a clinical whole-body scanner. In this work, PET stability of the GE SIGNA PET/MR was examined during simultaneous scanning of aggressive MR pulse sequences. Methods: PET performance tests were acquired with MR idle and during simultaneous MR pulsing. Recent system improvements mitigating RF interference and gain variation were used. A fast recovery fast spin echo MR sequence was selected for high RF power, and an echo planar imaging sequence was selected for its high heat-inducing gradients. Measurements were performed to determine PET stability under varying MR conditions using the following metrics: sensitivity, scatter fraction, contrast recovery, uniformity, count rate performance, and image quantitation. A final PET quantitative stability assessment for simultaneous PET scanning during functional MRI studies was performed with a spiral in-and-out gradient echo sequence. Results: Quantitation stability of a 68Ge flood phantom was demonstrated within 0.34%. Normalized sensitivity was stable during simultaneous scanning within 0.3%. Scatter fraction measured with a 68Ge line source in the scatter phantom was stable within the range of 40.4%-40.6%. Contrast recovery and uniformity were comparable for PET images acquired simultaneously with multiple MR conditions. Peak noise equivalent count rate was 224 kcps at an effective activity concentration of 18.6 kBq/mL, and the count rate curves and scatter fraction curve were consistent for the alternating MR pulsing states. A final test demonstrated quantitative stability during a spiral functional MRI sequence. Conclusion: PET stability metrics demonstrated that PET quantitation was not affected during simultaneous aggressive MRI. This stability enables demanding applications such as kinetic modeling.

Image-derived input function estimation on a TOF-enabled PET/MR for cerebral blood flow mapping.

15O-H2O PET imaging is an accurate method to measure cerebral blood flow (CBF) but it requires an arterial input function (AIF). Historically, image-derived AIF estimation suffers from low temporal resolution, spill-in, and spill-over problems. Here, we optimized tracer dose on a time-of-flight PET/MR according to the acquisition-specific noise-equivalent count rate curve. An optimized dose of 850 MBq of 15O-H2O was determined, which allowed sufficient counts to reconstruct a short time-frame PET angiogram (PETA) during the arterial phase. This PETA enabled the measurement of the extent of spill-over, while an MR angiogram was used to measure the true arterial volume for AIF estimation. A segment of the high cervical arteries outside the brain was chosen, where the measured spill-in effects were minimal. CBF studies were performed twice with separate [15O]-H2O injections in 10 healthy subjects, yielding values of 88 ± 16, 44 ± 9, and 58 ± 11 mL/min/100 g for gray matter, white matter, and whole brain, with intra-subject CBF differences of 5.0 ± 4.0%, 4.1 ± 3.3%, and 4.5 ± 3.7%, respectively. A third CBF measurement after the administration of 1 g of acetazolamide showed 35 ± 23%, 29 ± 20%, and 33 ± 22% increase in gray matter, white matter, and whole brain, respectively. Based on these findings, the proposed noninvasive AIF method provides robust CBF measurement with 15O-H2O PET.

NEMA NU 2-2012 performance studies for the SiPM-based ToF-PET component of the GE SIGNA PET/MR system.

PURPOSE: The GE SIGNA PET/MR is a new whole body integrated time-of-flight (ToF)-PET/MR scanner from GE Healthcare. The system is capable of simultaneous PET and MR image acquisition with sub-400 ps coincidence time resolution. Simultaneous PET/MR holds great potential as a method of interrogating molecular, functional, and anatomical parameters in clinical disease in one study. Despite the complementary imaging capabilities of PET and MRI, their respective hardware tends to be incompatible due to mutual interference. In this work, the GE SIGNA PET/MR is evaluated in terms of PET performance and the potential effects of interference from MRI operation. METHODS: The NEMA NU 2-2012 protocol was followed to measure PET performance parameters including spatial resolution, noise equivalent count rate, sensitivity, accuracy, and image quality. Each of these tests was performed both with the MR subsystem idle and with continuous MR pulsing for the duration of the PET data acquisition. Most measurements were repeated at three separate test sites where the system is installed. RESULTS: The scanner has achieved an average of 4.4, 4.1, and 5.3 mm full width at half maximum radial, tangential, and axial spatial resolutions, respectively, at 1 cm from the transaxial FOV center. The peak noise equivalent count rate (NECR) of 218 kcps and a scatter fraction of 43.6% are reached at an activity concentration of 17.8 kBq/ml. Sensitivity at the center position is 23.3 cps/kBq. The maximum relative slice count rate error below peak NECR was 3.3%, and the residual error from attenuation and scatter corrections was 3.6%. Continuous MR pulsing had either no effect or a minor effect on each measurement. CONCLUSIONS: Performance measurements of the ToF-PET whole body GE SIGNA PET/MR system indicate that it is a promising new simultaneous imaging platform.

Design Features and Mutual Compatibility Studies of the Time-of-Flight PET Capable GE SIGNA PET/MR System.

A recent entry into the rapidly evolving field of integrated PET/MR scanners is presented in this paper: a whole body hybrid PET/MR system (SIGNA PET/MR, GE Healthcare) capable of simultaneous acquisition of both time-of-flight (TOF) PET and high resolution MR data. The PET ring was integrated into an existing 3T MR system resulting in a (patient) bore opening of 60 cm diameter, with a 25 cm axial FOV. PET performance was evaluated both on the standalone PET ring and on the same detector integrated into the MR system, to assess the level of mutual interference between both subsystems. In both configurations we obtained detector performance data. PET detector performance was not significantly affected by integration into the MR system. The global energy resolution was within 2% (10.3% versus 10.5%), and the system coincidence time resolution showed a maximum change of < 3% (385 ps versus 394 ps) when measured outside MR and during simultaneous PET/MRI acquisitions, respectively. To evaluate PET image quality and resolution, the NEMA IQ phantom was acquired with MR idle and with MR active. Impact of PET on MR IQ was assessed by comparing SNR with PET acquisition on and off. B0 and B1 homogeneities were acquired before and after the integration of the PET ring inside the magnet. In vivo brain and whole body head-to-thighs data were acquired to demonstrate clinical image quality.